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via this fidelity function. We then show that the distances between quantum channels lead naturally 
to a new Fisher information which quantifies the ultimate precision limit in quantum metrology, 
the ultimate precision limit can thus be seen as a manifestation of the distances between quantum 
channels. We also show that the fidelity on quantum channels provides a unified framework for 
perfect quantum channel discrimination and quantum metrology, in particular we show that the 
minimum number of uses needed for perfect channel discrimination is exactly the counterpart of the 
precision limit in quantum metrology, and various useful lower bounds for the minimum number of 
uses needed for perfect channel discrimination can be obtained via this connection. 
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I. INTRODUCTION 

Fidelity, as a measure of the distinguishability between quantum states [THS], plays an important role in many areas 
of quantum information science, for example it is related to the precision limit in quantum metrology [ 4 ], serves as a 
measure of entanglement preservation through noisy quantum channels [ 5 ] , and a measure of entanglement preservation 
in quantum memory [ 6 ] ; it has also been used as a characterization method for quantum phase transitions [ 7 ] , and a 
criterion for successful transmission in formulating quantum channel capacities [ 8 ] . 

Unlike the fidelity of quantum states which is defined directly on quantum states, most commonly used measures 
for the distinguishability of quantum channels are defined indirectly through the effects of the channels on the states. 
For example the diamond norm, which is defined as \\Ki — i^2||o = ll^i ^ Ia{psa) — K2 0 7 a(ps'a)||i|SHII]( 

here ||X||i = TrVXtX, pSA denotes a state on system+ancilla, and I a denotes the identity operator on the ancillary 
system), is induced by the trace distance on quantum states \\pi — P2II1; another measure on quantum channels which 
is defined as arccosFmin(i^i, i^2) = arccosmiup^^ F.^lK^ a( pQ a), Ko^I a( pa A)] ^T^n^ , is induced by the fidelity on 

quantum states Fs{pi^p2) = Tr\Jpi p2pi . These induced measures through quantum states lack a direct connection 
to the properties of quantum channels, which severely restrict the insights that can be gained from these measures. 
A direct measure on quantum channels is expected to provide more insights thus highly desired. 

In this paper we provide a fidelity function defined directly on quantum channels, and show that this fidelity function 
on quantum channels, together with the classical fidelity on probability distribution and the fidelity on quantum states, 
form a hierarchy of fidelity functions in terms of optimization. This fidelity function on quantum channels also lead to 
various distance measures defined directly on quantum channels, in particular we show the Bures angle and the Bures 
distance can be extended to quantum channels. We then show the distance between quantum channels leads naturally 
to a new Fisher information on quantum channels which quantifies the ultimate precision limit in quantum metrology. 
We also show that this fidelity function provides a unified framework for perfect quantum channel discrimination and 
quantum metrology, in particular we show the minimum number of uses needed for perfect channel discrimination is 
exactly the counterpart of the precision limit in quantum metrology, and various useful lower bounds for the minimum 
number of uses needed for perfect channel discrimination can be obtained via this connection. 


II. FIDELITY FUNCTION ON QUANTUM CHANNELS 


We start by defining the fidelity function on unitary channels then extend it to noisy channels. 

For a m X m unitary matrix U, we denote as the eigenvalues of U, where Oj G (— tt, tt] for 1 < j < m and 

we call Oj the eigen-angles of U. We define(see also pTUT 6 ] i || U ||max= maxi<j<^ \ Oj |, and || U \\g as the minimum 
of II e'^^U Umax over equivalent unitary operators with different global phases, i.e., || U \\g= min^^]^ || e'^^U ||max- We 
then define 


C{U) = 


\\U\\g, if\\U\\g<^, 

f, if\\U\\A^- 


( 1 ) 


Quantitatively C{U) is equal to the maximal angle that U can rotate a state away from itself[T 6 l [171 El] 5 he., 
cos[C{U)] = min|^^ \ {2p\U\^)\. For mixed states it can be written as cos[( 7 (I/)] = miiip Fs{p^ UpU^). 

If ^max = 0 i > O2 > ' ’' > Om = ^min ^rc arranged in decreasing order, then C{U) = ^^ax-^min O^^ix — ^min < 

TT m- We then define &qc{Ui, U2) = C{U\U2)^ here Ui and U2 are unitary operators on the same Hilbert space(we 
can expand the space if they are not the same). It is easy to see that 


cos[Qqc{UuU2)]=cos[C{uIU2)] 

= rmnFs{U^pUlU2pUl), 

P 


( 2 ) 


©gc(1^15^/2) thus corresponds to the maximal angle between the output states of Ui and t/2(however we note that 
the definition of Qqc{Ui,U2) is independent of the states). We then denote Fqc{Ui,U2) = cos[ 0 gc'(I/i, I/2)] as the 
fidelity between Ui and I/2. For unitary channels this is equivalent to the fidelity function proposed previously in mi. 

We now generalize this to noisy quantum channels. A general quantum channel A, which maps from mi- to m2- 
dimensional Hilbert space, can be represented by Kraus operators, K{ps) = where 

Equivalently it can also be written as K{ps) = TrE{UEs{\^E){^E\ ^ Ps)U\^g)^ where |0£;) denotes some standard 
state of the environment, and Ues is a unitary operator acting on both system and environment, which we call as 
the unitary extension of K. 
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We define 0 gc'(^i,^ 2 ) = ^'^'^{UEsi^UES2}^QciUESi,UES2) and Fqc{Ki,K2) = cos 0 gc'(^i, ^2), where UESi 
are unitary extensions of i G { 1 , 2 }. In Appendix \A\ we show that the optimization can be taken by fixing one 
unitary extension and just optimizing over the other unitary extension, i.e., 

0Qc(^i,^2) = min &qc{UesiiUes2) 

( 3 ) 

= min &qc{UesiiUes2)- 

UeS 2 


In terms of Fqc{Ki, K2) it can be written as 

Fqc{Ki,K 2 ) = ^^'^Fqc{UeS 1 ,UeS 2 ) 

= maxFQc(C^£;5i Vbs2)- 

UeS 2 

This can be seen as the counterpart of Uhlmann’s purification theorem on quantum states [ 22 ] fhowever the proof 
does not use Uhlmann’s purification theorem [E]). In Appendix |B] we show that 0 gc'(i^i, i^2) is a metric and can 
be computed directly from the Kraus operators of Ki and K2 as [ 18 ] 


QQciKuK^) 


= arccos max -A^ 
ll^ll<i 2 


,{KwFKI)^ 


( 5 ) 


here Xmin{Kw + denotes the minimum eigenvalue of Kw + Kw with Kw = WijFl-F2j^ Fu and F 2 j denote 
the Kraus operators of Ki and K2 respectively, Wij denotes the ij-th entry of a q x q matrix W with ||IU|| < 1 where 
II • II is the operator norm which corresponds to the maximum singular value, here W arises from the non-uniqueness 
of the Kraus representations. Thus 


Fqc{Ki,K 2 ) 


max /'mm 
\\W \\<1 2 


^ '^min (Kw + kI). 


(6) 


We emphasize that Fqc is defined directly on quantum channels without referring to the states, such direct connection, 
in contrast to the induced measure, is crucial when applying the fidelity to channel discrimination and quantum 
metrology as we will show later. Furthermore the fidelity can be formulated as a semi-definite programming and 
computed efficiently as max||w||<i + F^w) ~ 


s.t. 


1 

max -t 

I VFt 
VF I 


^ 0 , 




w 


-tl hO- 


( 7 ) 


Analogous to the Bures distance on quantum states Bs{pi, P 2 ) = \/2 — 2 Fs{pi^ p 2 )^ we can similarly define a Bures 
distance on quantum channels as Bqc{Ki^ K2) = y /2 — 2 Fqc{Ki^ K2). In Appendix [ a| we prove an intriguing and 
useful connection between Bqc{Ki^ K2) and the minimum distances between the Kraus operators of Ki and K2 as 


BIc{Ki,K2)= ,min \\Y,{Fu - {Fu - F2i)\\ 
{Fii},{F 2 i] V 


where {Fu}, {F2i} are the sets of all equivalent Kraus representations of Ki and K2 respectively. This connection is 
particular useful in studying the scalings of the distance between quantum channels as we will show later. 

In which sense we call Fqc{Ki, K2) a fidelity function? It turns out that Fqc{Ki, K2) = min^^^ [i^i ( 8 ) 

Ia{psa)^ F2 0 Ia{psa)]‘ To see this, it is proved in the supplemental material of Ref. [ 18 ] that 

minF5[i^i 0 Ia{psa)^F 2 0 Ia{psa)] = max ^ '^min {Kw+Kl), (8) 

PSA \\W \\<1 Z 

which coincides with Eq. From this relationship it is also immediate clear that Fqc{Ki, K2) is stable, i.e., 

Fqc{Ki ^I,K2^I) = Fq(^Ki,K2). This result gives an operational meaning to Fqc{Ki, , K2). We emphasize that 
although we made connections between Fqc{Ki,K2) and the minimum fidelity of the output states, Fqc{Ki, K2) is 
defined directly on quantum channels and does not depend on the states. The definition and the operational meaning 
of Fqc{Ki, K2) play distinct roles in applications, the operational meaning provides a physical picture while the 
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direct definition brings insights which enable or ease the proofs and computations, which will be demonstrated in the 
applications. This is in analogy to how fidelity of quantum states is connected to the classical fidelity Fs{pi^p 2 ) = 
minj£;j Fc'(Pi,P 2 ): here Fc{pi^P 2 ) = \/^\/^ denotes the classical fidelity with pn = Tr{piEi) and p 2 i = 

Tr{p 2 Ei)^ {Ei} denotes a set of Positive Operator Valued Measurements(POVM) [3], here similarly the fidelity between 
quantum states has the operational meaning as the minimum classical fidelity, however the fidelity between quantum 
states is defined directly on quantum states which is independent of the measurements and such direct definition has 
provided numerous insights which would be hindered with just the classical fidelity. 

It is known that the trace distance and the fidelity between quantum states have the following relationships [19] 

1 - Fs{pi,p2) < - \\pi - P 2 II 1 < \/l “ E^{pi,p2), (9) 

from which it is straightforward to get the relationships between the diamond norm and the fidelity of quantum 
channels. This can be obtained by substituting pi = Ki 0 Ia{psa) and p 2 = K 2 ^Ia{psa)i then optimizing over psA 

maxi - Fs[Ki Ia{psa), ^2 ^ Ia{psa)] < max}-\\Ki ^ Ia{psa) -^2 ^ Ia{psa)\\i 

PSA PSA Z 

< max a/i - Eg[Ki (g) Ia{psa),K2 O Ia{psa)], 

PSA v 

which gives 

1-Fqc{K^,K2)<\\\K^-K2\U<^1-FIAKi,K2). ( 11 ) 

Since Fqc{Ki^K2) can be computed directly from the Kraus operators, this also provides a way to bound the diamond 
norm using the Kraus operators. 

In [20| the Choi matrices of the quantum channels are used to compute the fidelity between the channels, which 
corresponds to the fidelity between the output states of two quantum channels when the input state is taken as the 
maximal entangled state. As the maximal entangled state is in general not the optimal input state, the fidelity thus 
defined does not have operational meaning as the minimum fidelity of the output states, thus can not be related to 
the ultimate precision limit in quantum metrology etc (instead related to the precision limit when the probe state is 
taken as the maximally entangled state). 


III. A UNIFIED FRAMEWORK FOR QUANTUM METROLOGY AND PERFECT CHANNEL 

DISCRIMINATION 


Next we demonstrate the applications in quantum information science, in particular we show how the fidelity 
provides a unified platform for the ultimate precision in quantum metrology and the minimum number of uses needed 
for perfect channel discrimination. 

The task of quantum metrology, or quantum parameter estimation in general, is to estimate a parameter x encoded 
in some channel this can be achieved by preparing a quantum state psA and let it go through the extended channel 
Kx C) Ia with the output state px = Kx 0 Ia{psa)- By performing POVM, {Ey}, on px one gets the measurement 
result y with probability p{y\x) = Tr{Eypx). According to the Cramer-Rao bound [241127] . the standard deviation 
for any unbiased estimator of x is bounded below by Sx > , ^ =, where 6x is the standard deviation of the 

estimation of x, Jc[p{y\x)] is the classical Fisher information and n is the number of times that the procedure is 
repeated. The classical Fisher information can be further optimized over all POVMs, which gives 


Sx > 


1 

^nmaxiijj Jc\p{y\x)] 


1 

\/nJs{Pxy 


( 12 ) 


where the optimized value Js{px) is usually called the quantum Fisher information[4l l24l l25l [28] . here for distinguish 
we will call it the quantum state Fisher information. 

We first recall established connections between the fidelity functions and the Fisher information. Given px and 
its infinitesimal state px+dxi for a given POVM {Ey}^ the classical fidelity between p{y\x) = Tr{Eypx) and p{y\x + 
dx) = Tr{EyPx+dx) is given by Fc\p{y\x),p{y\x + dx)] = Y^y, \/p{yi\^)\/p{yi\^ + dx) which defines an angle as 
co^Qc[p{y\x)^p{y\x + dx)] = Fc[p{y\x)^p{y\x + dx)]. The classical Fisher information is related to the classical 
fidelity as \Jc[p{y\x)]dx‘^ = 2 — 2Fc[p{y\x)^p{y\x + dx)] up to the second order of dx[4], this can also be written as 


Jc\p{y\x)] 


A.el.\p{y\x),p{y\x F dx)\ 


(13) 
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If we optimize over {Ey} the classical fidelity then leads to the fidelity between quantum states as [4] 

F(j\TT{^EyPx\TT{^Eypx-\-dx)\ — Egi^Px-i Px-\-dx) 

\^ys 


and the classical Fisher information leads to the quantum state Fisher information Js{px) = ^^^{Ey} Jc[p{y\x)] and 
up to the second order of dxjH [28] 

~^Js{Px)dx =2 2Fs(^px’, Px-\-dx) - 

If we denote cos Ss{px^ Px+dx) = Fs{px, Px+dx)^ then 


^['^-COsSsipx^Px+dx)] 

d S\Px) — 7 9 

dx^O 

_ siPx’! Px+dx) 

dx^O dx^ 


(16) 


The precision can be further improved by optimizing over the probe states, which leads to the ultimate local 
precision limit of estimating x from K^. Intuitively, this ultimate precision limit should be quantified by the distance 
between Kx and its infinitesimal neighboring channel Kx+dx^ in a way analogous to how Bures distance of quantum 
states quantifies the precision limit of estimating x from the state px^- However although much progress has been 
made on calculating the ultimate precision limit [291137] . such a clear physical picture has still not been established 
after more than two decades since Braunstein and Caves’s seminal paper [4|, this is mainly due to the lack of proper 
tools on quantum channels. Here we show that the fidelity between quantum channels can be used to establish such 
a physical picture, which also leads naturally to a new Fisher information on quantum channel. 

Further optimizing over the probe states 

^^\js{px)dx^ = 2 -2xmnFs{px,Px+dx) 

PSA 4 PSA 

= 2-2FQc{Kx,Kx+dx) 

~ Fx+dx\ 

this leads naturally to a quantum channel Fisher information Jqc{Fx) = maxp^^ Js{px) which is similarly related to 
the distance on quantum channels as 


Jqc{K,) 


- lim 

dx^O 


^BlciKx 


Fx-\-dx^ 


dx‘^ 


lim Kx-\-dx)\ 

dx^Q dx‘^ 

T Fx-\-dx) 

hm —^- — -. 

dx^O dX"^ 


(18) 


The quantum channel Fisher information quantifies the ultimate precision limit upon the optimization over the 
measurements and probe states 


Sx > 


1 

\/nJQc{Kx) 


1 


Xn 


2QQc{Kx,KxJrdx) ’ 

\dx\ 


(19) 


This connects the precision limit directly to the distance between quantum channels which provides a clear physical 
picture for the ultimate precision limit. The scaling of the ultimate precision limit can now be seen as a manifestation 
of the scaling of the distances between quantum channels as we now show. 

Two schemes on multiple uses of quantum channels are usually considered in quantum parameter estimation, the 
parallel scheme and the sequential scheme as shown in Figj^ We will show that for both schemes, the scaling of the 
distances between two quantum channels are at most linear, which underlies the scaling for the Heisenberg limit. 

For parallel scheme with N uses of a channel K as shown in Fig|^ the total dynamics can be described by ^Ia> 

If we denote Ues as one unitary extension of then is a unitary extension of as shown in Fig|^ Given 

two channels Ki and K 2 , we choose Uesi and Ues2 as the unitary extension for Ki and K 2 respectively which satisfies 
Qqc{Ki^ K2) = &qc{UesijUes2)- Now as and unitary extensions of and respectively. 
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FIG. 1: (a) Parallel Scheme, (b)Sequential Scheme. 



FIG. 2: Parallel scheme with multiple 
uses of the channel. 


FIG. 3 : A unitary extension of the 
parallel scheme. 


we then have 

eQc{Kr,Kr) < eQc{u§^„u§g) 

= C[{Uis,UES2r^] ^ 20 ) 

< NCiUisiUES2) 

= NQqc{K^,K2). 

For the sequential scheme, we consider the general case that controls can be inserted between sequential uses of the 
channels. Any measurements that are used in the control can be substituted by controlled unitaries with ancillary 
systems, the controls interspersed between the channels can thus be taken as unitaries, which is shown in Figj^ 
Parallel scheme can be seen as a special case of the sequential scheme by choosing the controls as SWAP gates on the 
system and different ancillary systems [ 36 ]. We show that with N uses of the channel, the distance is still bounded 
above by NQqc{Ki,K2). 

We present the proof for the case of A" = 2 , same line of argument works for general N. For A = 2 , one 
unitary extension of U2K1U1K1 is U2Ue2SiUiUeiSi^ similarly U2UE2S2U1UE1S2 is a unitary extension of I/2A2I/1A2, 
here UEjSi denote a unitary extension of A^, i = 1,2, with Ej as the environment. We can choose UEjSi such 
that Qqc{Ki, K2) = QQc{UEjSi^UEjS2)^ here all operators are understood as defined on the whole space so the 
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FIG. 4: Sequential scheme with multiple uses of the channel. 



FIG. 5: A unitary extension of the sequential scheme. 


multiplication makes sense, for example the control I/i, which only acts on the system and ancillaries, is understood 
as Ui 0 Ie^ an operator on the whole space including the environment. We then have 

eQc{U2KiUiKi,U2K2UiK2) 

<Qqc{U2Ue2SiUiUeiSI,U2Ue2S2UiUeiS2) 

=C[Ui^SlUluisiUlU2UE2S2U,UE,S2] 

=C[U],^SlUlui^SlUE2S2UEE,S2] ( 21 ) 

<C[uIsiUe2S2]+C[uIsiUe,S2] 

=2eQciKuK2), 


i.e., with two uses of the channel, the distance is bounded above by 2&qc{Ki, K 2 ). With the same line of argument 
it is easy to show that with N uses of the channel the distance is bounded above by NQqc{Ki, K 2 ). 

Substitute Ki with and K 2 with Kx+dxi we have QQc{NKx^NKx+dx) < NQqc{Kx^ Kx+dx) for both schemes. 
From Eq.(19) the ultimate precision limit is then bounded by 


5x > 


1 


lim, 


dx^O 


2QQc(Kx,Kx+dx) 

\dx\ 




( 22 ) 


the scaling 1/N is called the Heisenberg scaling, which, as we showed, is just a manifestation of the fact that the 
distance between quantum channels can grow at most linearly with the number of channels. 

































































( 23 ) 


For N uses of the channels under the parallel scheme we can also obtain a tighter bound as 

2-2coseQc{Kf^,Kf^) 

< N\\2I -Kw-KI\\+ N{N - 1)||/ - Kwf, 

here Kw = ELi Ej=i F 2 j as previously defined, and the inequality holds for any W with ||W^|| < 1 (see 

Appendix [^. In the asymptotical limit, N{N — 1)||J — -AwlP is the dominating term, in that case we would like to 
choose a W minimizing ||/ — Kw\\ to get a tighter bound. This can be formulated as semi-definite programming with 
min|ny||<i \\I - Kw\\ = 


s.t. 


min 

t 

{I 

ITt 

V w 

/ 

( 

tl 

U- 

- Kw 


^0, 


(I 


- Kw)^ \ 
tl j 


0 . 


(24) 


If we let Ki = Kx and K 2 = then Eg. ([2^ provides bounds on the scalings in quantum parameter estimation, 

which is consistent with the studies in quantum metrology!^ [3Ql [32l l35l |36| but here with a more general context 
(see also Ref. [H]). 

Next we show how the tools unify quantum parameter estimation and the perfect quantum channel discrimination pTl 

EHHS]. 

Given two quantum channels Ki and they can be perfectly discriminated with one use of the channels if 
and only if there exists a psA such that Ki (g) Ia{psa) and K 2 G Ia{psa) are orthogonal, i.e., min^^^ (g) 

Ia{psa)^K2 G Ia{psa)] = 0, which is the same as 0qc(^i,^ 2) = f- When Ki and K2 can not be perfectly 
discriminated with one use of the channel, finite number of uses may able to achieve the task[42]. This is in contrast 
to the perfect discrimination of non-orthogonal states which always requires infinite number of copies. The minimum 
number of uses needed for perfect channel discrimination should satisfy &qc{NKi, NK 2 ) = f. The perfect channel 
discrimination is thus determined by the distances between quantum channels, and the scalings of &qc{NKi, NK 2 ) 
obtained before can be used to determine the minimum N. For example, from Qqc{NKi, NK 2 ) < NQqc{Ki, K 2 ) 
we can obtain a lower bound on N as 


^2eQc{K,,K2)^' 


(25) 


where [x] is the smallest integer not less than x. This bound is tighter than existing bounds for noisy channels}^ and 
for unitary channels it reduces to the formula which is known to be tight fTTI . For noisy channels under the parallel 
scheme we can also substitute QQc{Kf ^= f into the inequality (^23| to get a tighter bound. 

The lower bound on minimum N can also be obtained via a connection to quantum metrology. Given two channels 
Ki and i^ 2 , let K^^x G [a, b] as a path connecting Ki and K 2 . With N uses of the channel under the parallel strategy 

we have y Jqc{Kx^) = limda;_>o —~—^^e triangular inequality 




dx^o ax 


= l f ^JQc{Kf^)dx. 

^ J a 


(26) 


This connects the prefect channel discrimination to the ultimate precision limit. By choosing different paths various 
useful lower bounds on the minimum number of uses for perfect channel discrimination can be obtained. 

For example, given Kq{p) = and Ki = where cri,cr 2 and (73 are Pauli matrices and 

assume 0 = 0 . 8,77 = 0.5. For the parallel strategy the lower bound given by Eq.(25) is > \ 20 QciKo Ki) 1 ~ 
we choose a simple path Kx = (1 — x)Ko + xKi^ x G [ 0 , 1 ], which is a line segment connecting Kq to i^i, then with 
the connection provided by Eq.([^) we obtain N > 4. Other paths may be explored to further improve the bound. 
By using the inequality (23) with the W obtained from the semi-definite programming that minimizes ||/ — we 

get > 5. For any N we can also choose the W to minimize A^||2/ — Kw — + N{N — 1)\\I — it turns 

out the minimum N such that min||vj/||<i A^||2/ — Kw — Ky^W + N{N — 1)\\I — > 2 is 6 , thus N > 6. For 

comparison we also explicitly computed the actual distance GgciK^^^Kf^) with the increasing of A^, it turns out 
that the minimum N such that Gqc{K^^, Kf^) = ^ is actually 6 . All computations here are done with the GVX 
package in Matlab[44]. 
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IV. SUMMARY 

A fidelity function defined directly on quantum channels is provided, which leads to various distance measures 
defined directly on quantum channels, as well as a new Fisher information on quantum channel. This forms another 
hierarchy for fidelity functions and Fisher information as shown in the table: 

Fc{pi,P2) -^ Fs{pi,P2) -^ Fqc{Ki,K2) 

i Y Y 

0c(Pl,P2) -^ Qs{pl,p2) -^ Qqc{Ki,K2) 

Y Y Y 

Jc\p{y\'^y\ ^ Js{px) ^ Jqc{Kx) 

where cos©^ = Fi and Ji = hm^a,^o ^ {C, S', QC}. In this table the functions on quantum states equal 

to the optimized value over all measurements of the corresponding functions on probability distribution, and the 
functions on quantum channels equal to the optimized value over all probe states of the corresponding functions 
on quantum states. This framework connects quantitatively the ultimate precision limit and the distance between 
quantum channels, which provided a clear physical picture for the ultimate precision limit in quantum metrology. It 
also provide a unified framework for the continuous case in quantum parameter estimation and the discrete case in 
perfect quantum channel discrimination, with this framework the progress in one field can then be readily used to 
stimulate the progress of the other field. We expect these tools will find wide applications in many other fields of 
quantum information science. 


Appendix A: Formula to compute 0 qc(I^i,^2) 


We show that the distance between two quantum channels 0gc'(i^i, ^^ 2 ) = '^'^'^{UesiFes 2 }^Qc{Uesi^Ues 2 ) 
can be computed from the Kraus operators of Ki and K 2 as Qqc{Ki^K 2 ) = ^£;s' 2 ) = 

arccosmax|w|<i \^inin[Kw + ^ESi are unitary extensions of i G {1,2} and Xmin{Kw + ^w) denotes 

the minimum eigenvalue of Kw + Kw with Kw = ^ 2 j denotes the Kraus operators of Ki and 

K 2 , Wij denotes the ij-th entry of a g x g matrix W with ||IF|| < 1(|| • || is the operator norm which equals to the 
maximum singular value), g is the number of the Kraus operators. Furthermore the minimization on both Uesi and 
Ues 2 can be reduced to the minimization of just one 

Qqc{Ki,K 2 ) = min Qqc{Uesi,Ues 2 ) 

{UeS1,UeS2} 

= min 0qc{Uesi,Ues2) (A1) 

Uesi ^ ' 

= min Qqc{Ues\,^es 2 )‘ 

UeS2 


We start by a general unitary extension for any given channel K{p) 
maps from a mi- to m 2 - dimensional Hilbert space. 


Y!]=iPjpUj with E]=iFjP 3 = ’^'hich 


Ues = [We ( 8 ) Im2 


~Fi * * • • • 

P '2 * * • • • * 

) Fg * * 

0 * * 

_ 0 



u 


(A2) 


where We G U{p) only acts on the environment and can be chosen arbitrarily, here U{p) denotes the set of p x p 
unitary operators with p > g as p — g zero Kraus operators can be added. Here only the first mi columns of U are 
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fixed, the freedom of other columns can be represented as 


Ues = {We < 




Fi * * 
F2 * * 


* * 


0 **•••*. 


0 y 


(A3) 


where V can be any unitary. 

For two channels Ki and K 2 , with Ki{p) = ^ 2 {p) = unitary extensions 

be written as 


can 


Uesi = (Wei ^1712 ) 


'Fll * * 

Fi2 * * 

Fiq * * 

0 * * 


0 * * 


/mi 0 

0 Vi 


(A4) 


UeS2 = {We2 ^ /m2) 


'F 21 * * 
F 22 * * 

F2q * * 
0 * * 

_ 0 * * 


i-mi ^ 
0 ^2 


(A5) 


then 


UI;siUeS2 - 


ni 0 

_ 0 


K\y * * 
* * * 


* * * 

* * * 


Irm 0 

0 V 2 


here Kw = Eij FiEj, where Wij is the ij-th entry of IF, here W is the first q x q block of wI^We 2 , i-e., 

t [if *1 

WIj^We 2 = ^ II^11 — conversely for any IF with ||1F|| < 1 it can be imbedded as the 

first q X q block of a unitary matrix [45]. Thus by varying IFei and 1 Fe 2 we can take IF to be any q x q matrix with 
||fF|| < 1 . mmuEsi,UES 2 Qqc{Uesi,UeS 2 ) = C{Ul;g-^UES 2 ) is now reduced to the optimizing over Vi, 

V 2 and W. 

First note that for a fixed W, the first block of U^siUes 2 is always iFvv, as 


7mi 0 

0 y/ 


and 


‘-mi A/ 

0 ^2 


do not change 
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the first block. It has been shown in [ 46 ] that for any unitary that has Kw as the first block 


U 


■K^ 


VF * * 
* * * 


* * 
* * 


ll^llmax > arccos[^Amin(^vi/+^VF)], where 11 [/|| max is defined in Eq.(l) of the main text. Thus we have \\Uljg^UES2\\max > 
arccos[^Amin(^w + ^w)]' What’s more it was also shown that there exists a unitary V2 with 


Uv,= 


'K\y * * • • • *' 
>1< >1< . . . 


* * * 

* * * 


Im, 0 

0 ^2 


such that \\Uv2 Umax = arccos[|Aniin(-ft^VK + I^w)] achieves the bound [i6]. Similarly the bound can also be achieved by 
exploring the freedom in rows, i.e., there exists a unitary Vi, 


Uv, = 


/mi 0 

_ 0 vl_ 


’K\^ * * • • • * 


^ ^ >1< 


such that ||t/yj|max = arCCOs[|Amin(^lV+^Vv)]- Thus for a fixed Vk, minjy^^y^} Umax = arCCOs[^Amin(^VV + 


^w) 


Next we optimize over W. Basically we need to find W such that arccos ^Amin[^w + ^w] i^ minimized, which is 
equivalent to find max|w|<i |Amin[^w+^vv^]- fdat the freedom of global phase from || • ||max to || • ||^(see main text 
for definitions) has been included in the freedom of W and since max|w|<i ^Amin[^w+^vv] — ^^min[^o+^o] = 0? we 
have arccosmax||y|<i ^XminlKw+ ^ 1 ^] < f • Thus mmuEsi,UES2 C{Uljg^UES2) = arccosmax|^|<i ^Amin[^iv+ ^^^7], 


Qqc{Ki,K2) = min Oqc{Uesi,Ues2) 

UeS1,UeS2 


= arccos max -An 
\w\<i 2 




(A6) 




It is obvious that the freedom of W can be achieved by only varying Wi or IT2, thus the equality can be attained by 
just exploring the freedom of Vi and ITi, or V2 and IT2. We then have 


min Oqc{Uesi,Ues2) 

UeS1,UeS2 

= min &qc{UesiiUes2) 
Uesi 

= min &qc{UesiiUes2) 

UeS2 

= arccos max -Xmin\Kw + 
\W\<1 2 


(A 7 ) 
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Next we show that this distance measure has a connection to the minimum distance between equivalent Kraus 
operators. Given two quantum channels, Ki{ps) = ^uPs^l^ and K2{ps) = ^2iPsF2i{zeTO Kraus operators 

can be appended if the number of the Kraus operators are not the same), by appending additional p — q zero Kraus 
operators, we have the Kraus operators for Ki and K2 as {^11,^12, * • * , 0 , • • • , 0 } and {F21, ^22, * * * , F2q^ 0 , • • • , 0 } 

respectively. Equivalent Kraus operators for Ki and K2 can be represented as Fu = and F2i = '^i/c^2/c 

where Uik and Vik are entries of [/, K e U(p) respectively, here I < i < p. Then 


.min \\Y{Fu - F2i)\Fu - F2i)\\ 
= .min \\ 2 I-j 2 iFlAi+F^,Fu\\ 


(A 8 ) 


where Kw = THj WijFl^F2j and Wij is the ij-th entry of IT, which is the first q x q block of U^V and can be any 
q X q matrix with || 1 T|| < 1 by varying U and T, i.e., by varying the equivalent representations of Ki and K2. Thus 


.min. \\y^{Fu-F2p{Fu-F2i)\\ 

(A 9 ) 

=2 - I max^ XminiKw + ATJv) 

=2 - 2 cos ©^^(Ari, ^2), 


we then have 


p 

BIc{K^,K2) = 2 - 2 cosQqc{K^,K 2 )= min \\T{Fu - F2i)HFu - F2i)\\. (AlO) 


Appendix B: ©qc(Ki,K 2) defines a metric on quantum channels 

We show that Qqc{Ki^K 2 ) defines a metric on quantum channels. 

First we show that 0Qc(t^i, U 2 ) = C{U\U 2 )^ where C is defined in the main text, is a metric on unitary channels. 
We start by listing some useful properties of C{U): 

C{V^UV) = C{U)-, 

C{Ui ® U 2 ) < C{Ui) + C{U 2 ); (Bl) 

C{UiU2) < C{Ui) + C{U2)\ 

where V is any unitary operator. The first equality is obvious from the definition; the second inequality can be easily 
verified using the formula 0(17) = when 6>niax —^min < f, the equality is saturated when 0 ( 1 / 1 ) 70 ( 1 / 2 ) < f; 

proof of the third inequality can be found in gaiiH]. 

It is obvious that Qqc{U, U) = 0 and 0Qc(h^i V 2 ) = 0 Qc(C^ 2 , Ui) > 0 if {7i 7 ^ U 2 - And since 

OQciUuUs) = CiUlUs) 

= CiUlU2U7s) (J32) 

<c(c/|c/ 2 ) + c(c/]c/ 3 ) 

= eQc(t/l V 2 ) + 0Qc(C/2 V 3 ). 

where for the inequality we have used the property that 0(UiU2) < 0(Ui) 7 0 ( 1 / 2 ). This shows that Oqc(Ui, U 2 ) is 
a metric on unitary operators. 

For two general channels, Qqc{Ki,K 2 ) = ToamuEsi^Qc{UESi,UES 2 ) = ToamuES 2 ^Qc{Uesi,UeS 2 ) where Uesi 
and Ues 2 are unitary extensions for Ki and K 2 respectively. It is easy to see that Qqc{Ki, K 2 ) = Qqc{K 2 , Ki) > 0 
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and the equality is saturated only when Ki = K 2 . We show that 0gc also satisfies the triangular inequality as 


eQc{K,,Ks) 

= min Oqc{Uesi,Ues3) 

Uesi 

= min C{Ul;g^UES3) 

Uesi 

= min C{Ul,g-JJES2Ul:g2UES3) 

Uesi 

< min [C{Uig^UES 2 ) + C{U]gg2UES3)] 

u ESI 


(B3) 


the last inequality is valid for any Ues 2 , specially we can choose the Ues 2 which minimizes C{U\^g 2 UES 3 ), thus 

Qqc{K^E3) < min C{uIs^Ues 2 ) + QQciK^Ei) 

Uesi (B4) 

= Qqc{Ki,K 2) + eQc{K2,K3)- 

Qqc{Ki,K 2 ) thus defines a metric on the space of quantum channels. 


Appendix C: Upper bound of the distance with N parallel channels 


Given two quantum channels, Ki{ps) = ^upsP^- and K2{ps) = appendingp—g zero Kraus 

operators, we have the Kraus operators for Ki and K 2 as {Fii, F 12 , * * * , 0, • • • ,0} and {F 21 , ^ 22 , * * * , ^ 2 g, 0, • • • ,0} 

respectively. All the equivalent Kraus operators for Ki and K 2 can be represented as Fu = UikFik and F 2 i = 
'^ikF 2 k where Uik and Vik are entries of U, K G U(p) respectively, here I < i < p. 

With N channels in parallel, one representation of the Kraus operators for can be written as - = 

Pii! ®Pii 2 similarly for Kf^ we have F2i^,2h,--- , 2 i« = Fi! F! = ELi -^iikEk 

are one Kraus operators of the l-th channel of similarly F^^ = '^iz/c^ 2 /c are one Kraus operators of the 

l-th channel of As -d^iv} and {F 2 i^^ 2 i 2 ,--- ,2iN} are just one particular Kraus representation of 

and respectively, we then have 


2-2cos0gc(i^f 


0Ar 


K^n 


)<ll E Fn, 


li2?"' d^AT 


-F, 


2ii ,2*2 J • ■ ■ d* Af 


}HFie 


li2j"' d^Af 


-F, 


2ii ,2i2 ,• • • ,2iiv / || i 


Since 


Fu 


lii,li2,--- ,1*A 


Foi 


=F^! 5 

di) 


— ^2ii 


2ii ,2^2, • • • ,2iE 
> 


^2i2 


I 

1*2 


Hn 2ii 

I r(i) 


^2iN 




by induction it is then easy to get that 

Fi 


odd?"' d*Af 






db 


2ii ,2^2, • • • ,2iE 


F, 

(2) 


Hn 


( 3 ) 






’ ^2i2 


(X-dd) 


l*Ar 


^2i2 


idW 

l*Ar 

,d^) 


(Cl) 


)d^^i 

' -^2iN J’ 


(C2) 


+dd 


2^2 


0 


1*7V 


p(N). 

^2iN ) 
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Thus 


1=1 


^ ^ jl^AT ~ -^2ii,2i2,-" 52i7v ) 4*iv ~ -^2ii,2i2,-" j2iAr) 

il,22,-" 4iV 
N 

=/ 0 j 0 • • • 0 

1 ii 

N h-l 

E E ® ® ® ® [E 

/2 = lZi = l ii-^ ii2 

I 

= ^/ 0j0---0{2J- iiTvi/ - 0 / 0 • • • 0 / 

• 0 - /) 0 ■ • ■ 0 (/ - 0 • ■ • 0 j + /i.c]. 


2^2 


AT 


Z=1 


) I + h.c] 


(C3) 


N h-l 

■E Ei^®- 

l2 = lh = l 

here again W is the first q x q block of WV and Kw = Yli=i Si=i with Wij as the ij-th entry of W. We 


then have 


2-2coseQc{Kf^,Kf^) 

— II ^ ^ (-^lii,li2,-" ,l*iv ~ -^2ii,2i2,-" ,2i7v)^(-^ln,l*2,---,1*JV ~ -^2ii,2i2,"-, 


2iAr j 


(C4) 


21,22,--- ,^7V 


< N\\2I -Kw-KI\\+ N{N - 1)IIJ -Kwf. 


[1] C. A. Fuchs and C. M. Caves, Phys. Rev. Lett. 73 , 3047 (1994), URL http://link.aps.org/doi/10.1103/PhysRevLett. 
73.3047 

[2] R. Jozsa, Journal of Modern Optics 41, 2315 (1994), http://dx.doi.org/10.1080/09500349414552171, URL http://dx. 
doi.org/10.1080/09500349414552171 

[3] C. A. Fuchs, arXiv pp. quant-ph/9601020 (1996). 

[4] S. L. Braunstein and C. M. Caves, Phys. Rev. Lett. 72, 3439 (1994), URL http://link.aps.org/doi/10.1103/ 
PhysRevLett.72.3439 

[5] B. Schumacher, Phys. Rev. A 54 , 2614 (1996), URL http://link.aps.org/doi/10.1103/PhysRevA.54.2614 

[6] K. Surmacz, J. Nunn, F. C. Waldermann, Z. Wang, 1. A. Walmsley, and D. Jaksch, Phys. Rev. A 74 , 050302 (2006), URL 
http://link.aps.org/doi/10.1103/PhysRevA.74.050302 

[7] S.-J. Gu, Int. J. Mod. Phys. B 24, 4371 (2010). 

[8] H. Barnum, E. Knill, and M. A. Nielsen, IEEE Trans. Info. Theor. 46, 1317 (2000), quant-ph/9809010. 

[9] A. Kitaev, Russian Mathematical Surveys 52 , 1191 (1997), URL http://www.turpion.org/php/paper.phtml7journal_ 
id=rm&paper_id=2155 

[10] A. Y. Kitaev, A. H. Shen, and M. N. Vyalyi, Classical and Quantum Computation (American Mathematical Society, 
Boston, MA, USA, 2002), ISBN 0821832298. 

[11] J. Watrous, arXiv p. 0901.4709 (2009). 

[12] A. Gilchrist, N. K. Langford, and M. A. Nielsen, Phys. Rev. A 71 , 062310 (2005), URL http://link.aps.org/doi/10. 
1103/PhysRevA.71.062310 

[13] V. P. Belavkin, G. M. DAriano, and M. Raginsky, J. Mathematical Physics 46 , 062106 (2005). 

[14] H. E. Chau, Quant. Inf. Comput. 11, 0721 (2011). 

[15] C.-H. E. Eung and H. E. Chau, Phys. Rev. A 88, 012307 (2013), URL http://link.aps.org/doi/10.1103/PhysRevA.88. 
012307 

[16] C.-H. E. Eung and H. E. Chau, Phys. Rev. A 90 , 022333 (2014), URL http://link.aps.org/doi/10.1103/PhysRevA.90. 
022333 

[17] A. Acin, Phys. Rev. Lett. 87 , 177901 (2001), URL http://link.aps.org/doi/10.1103/PhysRevLett.87.177901 

[18] H. Yuan and C.-H. E. Eung, to appear in npj Quantum Information (2017). 

[19] C. A. Euchs and J. van de Graaf, IEEE Trans. Inf. Theory 45, 1216 (1999). 

[20] Maxim Raginsky, Physics Letters A, 290 , 11-18 (2001). 

[21] G. M. D’Ariano, P. Lo Presti, and M. G. A. Paris, Phys. Rev. Lett. 87 , 270404 (2001), URL http://link.aps.org/doi/ 
10.1103/PhysRevLett.87.270404 

[22] A. Uhlmann, Rep. Math. Phys. 9 , 273-279 (1976). 

[23] H. Yuan and C.-H. E. Eung, arXiv p. 1506.01909 (2015). 

[24] C. W. Helstrom, Quantum Detection and Estimation Theory (Academic Press, 1976). 

[25] A. S. Holevo, Prohahilistic and Statistical Aspect of Quantum Theory (North-Holland, 1982). 

[26] H. Cramer, Mathematical Methods of Statistics (Princeton University,, 1946). 



15 


[27] C. R. B. Rao, Calcutta Math. Soc. 37 , 81 (1945). 

[28] S. L. Braunstein, C. M. Caves, and G. J. Milburn, Annals of Physics 247 , 135 (1996). 

[29] A. Fujiwara and H. Imai, Journal of Physics A: Mathematical and Theoretical 41, 255304 (2008). 

[30] B. M. Escher, R. L. de Matos Filho, and L. Davidovich, Nature Physics 7 , 406 (2011). 

[31] M. Tsang, New J. Phys. 15, 073005 (2013). 

[32] R. Demkowicz-Dobrzanski, J. Kolodynski, and M. Guta, Nature Communications 3 , 1063 (2012). 

[33] S. Knysh, V. N. Smelyanskiy, and G. A. Durkin, Phys. Rev. A 83, 021804 (2011), URL http://link.aps.org/doi/10. 
1103/PhysRevA.83.021804 

[34] S. Knysh, E. Chen, and G. Durkin, arXiv p. 1402.0495 (2014), URL http://arxiv.org/abs/1402.0495 

[35] J. Kolodynski and R. Demkowicz-Dobrzanski, New J. Phys. 15, 073043 (2013). 

[36] R. Demkowicz-Dobrzanski and L. Maccone, Phys. Rev. Lett. 113 , 250801 (2014), URL http://link.aps.org/doi/10. 
1103/PhysRevLett.113.250801 

[37] S. Alipour, M. Mehboudi, and A. T. Rezakhani, Phys. Rev. Lett. 112, 120405 (2014), URL http://link.aps.org/doi/ 
10.1103/PhysRevLett.112.120405 

[38] R. Duan, Y. Feng, and M. Ying, Phys. Rev. Lett. 98, 100503 (2007), URL http://link.aps.org/doi/10.1103/ 
PhysRevLett.98.100503 

[39] R. Duan, Y. Feng, and M. Ying, Phys. Rev. Lett. 100 , 020503 (2008), URL http://link.aps.org/doi/10.1103/ 
PhysRevLett.100.020503 

[40] C. Lu, J. Chen, and R. Duan, Quantum Info. Comput. 12, 138 (2012), ISSN 1533-7146, URL http://dl.acm.org/ 
citation.cfm?id=2231036.2231045 

[41] G. Chiribella, G. M. D’Ariano, and P. Perinotti, Phys. Rev. Lett. 101, 180501 (2008), URL http://link.aps.org/doi/ 
10.1103/PhysRevLett.101.180501 

[42] R. Duan, Y. Feng, and M. Ying, Phys. Rev. Lett. 103 , 210501 (2009), URL http://link.aps.org/doi/10.1103/ 
PhysRevLett.103.210501 

[43] A. W. Harrow, A. Hassidim, D. W. Leung, and J. Watrous, Phys. Rev. A 81 , 032339 (2010), URL http://link.aps.org/ 
doi/10.1103/PhysRevA.81.032339 

[44] M. Grant and S. Boyd (2011), URL http://cvxr.com/cvx/ 

[45] M.D. Choi, and C.K. Li, J. Operator Theory. 46, 435 (2001). 

[46] C.-H. F. Fung, H. F. Chau, C.K. Li, and N.S. Sze, Quantum Information and Computation 15, 0685-0693 (2015). 

[47] A. Childs, J. Preskill, and J. Renes. Quantum information and precision measurement. Journal of Modern Optics, 47(2- 
3):155-176, 2000. 

[48] H. Yuan, C.-H. F. Fung, Phys. Rev. Lett. 115, 110401 (2015). 



